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Abstract 

The harmonic polylogarithms (hpl's) are introduced. They are a generalization of Nielsen's 
polylogarithms, satisfying a product algebra (the product of two hpl's is in turn a combination 
of hpl's) and forming a set closed under the transformation of the arguments x — l/z and 
X = (1 — + The coefficients of their expansions and their Mellin transforms are harmonic 
sums. 
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1 Introduction. 

Euler's dilogarithm appeared very soon, even if with a different name, in the evaluation of radiative 
corrections in QED. The first occurrence is perhaps in the 1934 paper by G. Racah on the radiation 
by fast particles Q, whose function F(x) is equal to — Li2(— x) in Euler's notation. Two loop 
calculations |2| required the polylogarithms, Nielsen's generalization of Euler's dilogarithm. 
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More bibliographical indications as well as many relevant results are contained in the popular book 
by Lewin (note the change in the titles of the two editions of the book) . 

While the polylogarithms are the natural analytical tool to use when dealing with the (relatively) 
simple integrals appearing in calculations with a few loops, it is known that they will not be 
sufficient when the number of loops will be larger than has been considered thus far or when 
several different scales are present. In a recent publication the set of polylogarithms has been 
extended into something called 'multidimensional polylogarithms' [^. These functions seem to 
be very useful when more than one dimensionful parameter is involved. In principle they are a 
direct generalization of the definition of the power-series expansion of the polylogarithms to a 
multiparameter space. 

Besides the dilogarithm, Euler studied also harmonic sums. A recent publication by one of 
us |6| investigated harmonic sums and their applicability, in particular to formulas in Mellin space. 
These harmonic sums seem to be the natural functions for the results of moment calculations of 
deep inelastic structure functions when only massless quarks are involved]^. If indeed all these 
moments can be expressed in terms of harmonic sums, the class of functions that will represent the 
results in the regular x-space will be formed by the inverse Mellin transforms of these harmonic 
sums. In ref Q it was indicated how one could obtain at least numerical representations of these 
functions by means of numerical integration. 

In the current paper we study these functions in a more systematic way. We start with a 
recursive integral definition of a class of functions, which we will call the harmonic polylogarithms 
(hpl's), which are by construction a generalization of Nielsen's polylogarithms; it turns out, further, 
that an important subset of the hpl's is also a subset of the multidimensional polylogarithms 
of ref Q. Then we will study a number of their properties, including expressions for products 
of harmonic polylogarithms with the same argument, the behaviour at x = 0, 1, the relevant 
expansions around those points, the algebra of the hpl's and the identities between hpl's of related 
arguments. Then we study special values and numerical evaluation. Finally we study the Mellin 
transforms of the harmonic polylogarithms and find that indeed they give the harmonic sums and 
that there is a one to one correspondence between them. As a consequence the investigation also 
leads to a rather simple algorithm for the inverse Mellin transform, even though in general the 
length of the resulting formulae requires a computer implementation for dealing with the great 
number of terms which are generated. 

All algorithms that we present have been programmed in the language of FORM 0]. The 
resulting procedures can be obtained from the second author. 



^This can be shown for all two loop calculations to any order in the expansion parameter e. For three loop 
calculations such results do not exist yet, but a recent result by Broadhurst and Kreimer shows that only at the 
7-loop level the counter terms in the QCD beta function contains non-zeta like constants. 
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2 Definitions. 



The harmonic polylogarithms of weight w and argument x are identified by a set of w indices, 
grouped into a tD-dimensional vector and are indicated by H(m^(,; x). 
More expHcitly, for w = 1 one defines 

H(0; x) = Inx , 

H(l;x) = / ^ = -ln(l-x), 

Jo I - x' 

H(-l;x) = r^=ln(l+x). (1) 
For their derivatives, one has 

^R{a;x) = f{a;x) , (2) 
where the index a can take the 3 values 0, +1, —1 and the 3 rational fractions /(a; x) are given by 

/(0;x) = -, 
x 



-I ) 

1 — x 



/(-i;^) = (3) 

Note the (minor) asymmetry of Eq.(|^), in contrast with the higher symmetry of Eq.(|^). 

For w > 1, let us elaborate slightly the notation for the u;-dimensional vectors m^. Quite in 
general, let us write 

rhw = (a,m^-_i) , (4) 

where a = is the leftmost index (taking of course one of the three values 0, 1, —1), and rhyj-i 
stands for the vector of the remaining {w — 1) components. Further, 0^, will be the vector whose 
w components are all equal to the index 0. The harmonic polylogarithms of weight w are then 
defined as follows: ^ 

w 

while, if rhw ^ 0, 



H(0^;x) = — In^-x , (5) 



H(m^;x) = / dx' f{a;x') H(m^_i;x') . (6) 
Jo 

Quite in general the derivatives can be written in the compact form 

^H(m^; x) = f{a; x)H(m^_i; x) , (7) 
where, again, a = is the leftmost component of niw 
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In analogy with Eq.(|5|), if 1^, (— 1)^„ are the vectors whose components are ah equal to 1 or — 1, 
we have by applying recursively the definitions 



H(l^;x) = — (-ln(l-x)r , 
w] 

H((-lL;x) = ^ln-(l + x) . 



(8) 



Let us now have a look at the first few values of the indices. For w = 2 one has the 9 functions 



H(0,0;x) 
H(0,l;x) 
H(0,-l;x) 
H(l,0;x) 
H(l,l;x) 
H(l,-l;x) 
H(-l,0;x) 
H(-l,l;x) 
H(-l,-l;x) 



lln^x 



r^H(l;x') = - r^ln(l-x') 
JO x' Jo x' 

r^H(-l;xO= r^Hi + x') 

Jo x' Jo x' 



^ dx' 





^ dx' 



In x' 



Jo I — x' Jo I — x' 

/ f^H(-l;x')= / -^ln(l+x') 
Jo ^ - x' Jo 1 — x' 



^ dx' 
1 - s 
r dx' 

rx i-x fj^i 



^ dx' 

V- 

" dx' 
1 + 
^ dx' 



In x' 



(9) 



Those 9 functions can all be expressed in terms of logarithmic and dilogarithmic functions; indeed, 
if 

'■^ dx' 



j-X A^' 

Li2(x) =- / -^ln(l-x') 
Jo x' 

is the usual Euler's dilogarithm, one finds 



(10) 



H(0, l;x 
H(0,-l;x 
H(l,0;x 

H(l,l;x 

H(l,-l;x 
H(-l,0;x 
H(-l,l;x 

H(-l,-l;x 



Li2(x) , 
-Li2(-x) , 

— lnxln(l — x) + Li2(x) , 



^ln2(l-x) , 



Li2 



1 — X 



lnxln(l + x) + Li2(— x) , 
iln2(l + x) . 



In21n(l - x) - Li2 
Li2(-x) , 

ln21n(l + x) - Li2 



(11) 
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Something similar happens for harmonic polylogarithms and Nielsen's poly logarithms of weight 3; 
that is no longer true however from weight 4 on. To make an example, 

H(-l,0,0,l;x) = / -— ^Li3(x') (12) 
Jo 1 + x' 

cannot be expressed in terms of Nielsen's polylogarithms of the same weight, even allowing for 
slightly more general arguments (i.e. when considering, besides x, also — x, (1 + a;)/2, (1 — x)/2 
etc.). In other words, the set of the S"' harmonic polylogarithms of weight w is in general a much 
wider set of functions than the set of the Nielsen's polylogarithms. 

It follows from the definition that if ^ 0^ the hpl's vanish at x = 0: 

H(m^;0)=0, fh^^Qu)- (13) 

Likewise, if the leftmost index niyj is not equal to 1, {myj ^ 1), H(m^„; 1) is finite; it is also finite 
when fhuj = 1, but all the remaining indices fh^-i are zero, (mu,_i = 0^_i). In the remaining 
cases, i.e. niw = 1 and m^_i ^ Ow-i, II(m^;x) has a logarithmic behaviour at x = 1: more 
exactly, if the p leftmost indices are all equal to 1, II(m^; x) behaves for x ^ 1 as a combination of 
powers of ln(l — x) ranging from the maximum value p down to (the maximum power is decreased 
to p — 1 if the remaining w — p indices are all equal to zero; the study of the detailed logarithmic 
behaviours at x = 0, 1 will be carried out in Section 3). 

In dealing with specific cases and except for the smallest values of w, specifying explicitly all the 
components of m becomes quite cumbersome, so that a more compact notation is welcome. In the 
case that we ignore the functions of which the last index is zero we can use the same compactified 
notation as in ref This is to say that, proceeding from right to left, all zeroes are simply 
eliminated by adding at the same time one to the absolute value of the previous index to the right, 
as in 

H(0,0,l,0,-l;x) =H3,_2(x). (14) 

In terms of this notation and excluding, as already stated, the cases in which the rightmost index 
is zero, one can formulate the following: 

theorem: If mi ^ one has 

Hm,,.,mi(-X) = i-l)PR_m„-,-m,{x) ■ (15) 

The proof goes by induction and follows rather trivially from the definition of H. In the case that 
we use the notation in which the nii only have the value 0, 1, —1 the power of —1 is the number of 
indices that are not zero. 

In general we will write the indices of the H-functions as subscripts when we use the notation of the 
r.h.s. of Eq.(ll), while we will use the notation of the l.h.s. when the indices are supposed to have 
the values 0, 1, —1 only. In that last notation, to see the relation with the polylogarithms of Nielsen 
Sn,p{x), defined in let us indicate with 0„,lp as usual, two n-dimensional and p-dimensional 
vectors whose components are all equal to and 1 respectively; one then has 

Sn,p{x) = II(On, Ip] x) . (16) 
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As an obvious extension of the terminology, the product of two H- functions of weight wi and W2 
will be said to have total weight w = 'Wi + W2- In the following we will often encounter homogeneous 
"identities of weight w" , i.e. relations (or identities) involving the sum of several terms, where each 
term is equal to the product of an integer or rational fraction times a H-function of weight w or a 
product of several H-functions separately of lower weight but with total weight w. 

While the H-functions of weight w are linear independent, the same is not true for the wider 
set of all the homogeneous expressions of weight w. The redundance can be used for establishing 
a number of (homogeneous) identities expressing a H-function of some argument and weight w 
as a homogeneous expression of the same weight involving H-functions of the same or of related 
arguments (including constant arguments, such as for instance -|-1 or —1). The identities can be 
useful, typically, for exhibiting explicitly the behaviour at particular points (such as the logarithmic 
behaviour at or ±1) or for obtaining relations between H-functions of special arguments. Quite in 
general, while establishing such identities can be more or less wearisome, there is almost always a 
straightforward "standard method" for checking a given identity: one first verifies that the identity 
holds for a particularly convenient choice of the variable (or variables) and then differentiate it 
with respect to one of the arguments. In so doing one obtains another relation, however of lower 
weight, according to Eq.(^); the procedure can be iterated until a relation of weight 1 is eventually 
obtained, whose check is trivial (because the H-functions of weight 1 are just logarithms). 

Likewise, also the mathematical constants corresponding to the particular values of the H- 
function of weight w (such as the values at a; = 1 when finite) can be given the same weight w. 
Those values, aX x = ±1 or other simple arguments, are of particular interest by themselves, as it 
turns out that they can be expressed in terms of a very small number of mathematical constants, 
such as Riemann (^-functions. In 2 etc. We will see that they are connected to the sums to infinity of 
ref which have been systematically evaluated and tabulated^by one of us (J-V.) to weight = 9 
and can be evaluated basically to any weight, given enough computer resources^. In similar ways 
these sums have been evaluated under the name of Euler/Zagier sums by the authors of ref 
Hence, whenever H-functions aX, x = 1 will appear in this paper they can be regarded as known 
from ref Q or ref Q, provided their weight is not too large. It will also be shown that one may 
alternatively consider them as unknown constants, to be expressed in terms of that much smaller 
number of mathematical constants by systematically exploiting the many identities among H's of 
various arguments established in the rest of this paper. 



^ in ref [g] this was done only to weight = 7 

^An ahernative method to obtain the finite constants consists of their numerical evaluation to high precision and 
then fitting them to a presumed basis. Using this method Broadhurst |ll| has evaluated all finite objects at weight 
= 9 and some objects at the weights 10 and 11 
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3 Identities between functions of the same argument. 



Let us start by the integration by parts (ibp) identities. From the very definition, 

H(mi • • • mg; x) = / dx' /(mi; x')H(m-2 • • • m^; x') 
Jo 

rx 

= H(mi; x)H(m2 • • • m„; x) — / dx' H(mi; x')/(m2; x')H(m3 • • • m„; x') 

Jo 

= H(mi; x)H(m2 • • • m^; x) — H(m2mi; x)H(m3 • • • m^; x) 

+ H(m3m2m,i; x)H(m4 • • • m^; x) — • • • — (— l)^H(mg • • • mi; x) . (17) 

The above identity can be immediately verified, independently of its derivation, by the 'standard 
methods': it holds at x = 0; when differentiating with respect to x, one obtains a number of terms 
which are immediately seen to cancel out pairwise; therefore, the relation is true. This relation 
shows that in the case that rhq is symmetric and q is even the H-function reduces to products of 
lower weight functions. In general the relation can be used when it is important to reduce the 
number of H-functions with the highest weight as much as possible. 

Another important set of identities expresses the product of any two H-functions of weight wi 
and W2 as a linear combination of H-functions of weight w = wi + W2- Let us start from the case 
wi = 1; the identity reads 

H(a; x)H(mp, • • • , mi; x) = H(a, m^ • • • , mi; x) 

+ H(mp, a, mp_i • • • , mi; x) 

+ H(mp,mp_i,a,mp_2 • • •mi;x) 

+ ••• 

+ H(mp, • • • ,mi, a;x) . (18) 

It can be established by induction in p. For p = 1 it is almost trivial, corresponding to Eq.([T7|) 
for q = 2. Assume then that it holds for p — 1; take the identity for p — 1, multiply by /(mp; x) 
and integrate over x. In the r.h.s. we can do the integral and obtain all necessary terms except for 
the one starting with a. The l.h.s. can be integrated by parts to give the proper l.h.s. term plus 
another term that can be integrated and gives indeed the missing term. This completes the proof. 

Again, once established the identity can also be verified by the 'standard method': it holds at 
X = 0; the x-derivative consists of two groups of terms, a first group with the coefficient /(a;x) 
contains just two terms which cancel out immediately, plus a second group proportional to /{nip] x), 
which is nothing but the same relation at level p—1, so that the procedure can be repeated p times 
until everything cancels out. 

There is only one complication with Eq.(^). This concerns points in which one of the objects 
involved is divergent. Hence one cannot apply this equation for x = 1 in the case that either a = 1 
or nip = 1. This is explained better in the section on the algebraic properties. 

Eq.(|l8|) can be generalized to the product of two H-functions H(p; x)H(^; x); if p,q are the 
dimensions of p, q (or, which is the same, the weights of the two H-functions), the product is equal 
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to the sum of {p + qY-fplq^- terms, each term being an H-function of weight {p + q) with coefficient 
+1, obtained by choosing p indices in all possible ways (hence the binomial coefficients) and filling 
them from left to right with the components of p without changing their order, while the remaining 
q places contain the components of q, again without altering their order. This can be expressed 
with the formula 

R{p;x)R{q;x) = ^ R{f;x) (19) 

r=jAHq 

in which p^S q represents all mergers of p and q in which the relative orders of the elements of p 
and q are preserved. 

As an example, for p = 2,p= (a, b) and q = 3,q= {r, s, t) one has 

tl{a,b; x)}l{r, s,t; x) = }l{a,b,r,s,t;x)+}i{a,r,b,s,t;x) 

+ }l{a,r,s,b,t;x)+}i{a,r,s,t,b;x) 

+ H(r, a, 6, s, i; x) + H(r, a, s, 6, t; x) 

+ H(r, s, a, 6, t; x) + H(r, a, s, i, 6; x) 

+ H(r, s, o, t, 6; x) + H(r, s, t, a, 6; x) , (20) 

as can be easily checked, again, by the 'standard method'. 

The product identities Eq.(|l9|) can be used to single out the terms in ln(x) from H- functions 
whose indices have trailing (or rightmost) indices equal to zero (as we will see in the next section 
H-functions with no trailing zeroes can be expanded in series of x around x = 0, while H-functions 
with trailing zeroes develop logarithmic singularities at that point). For a = in Eq.([l^), recalling 
H(0;x) = ln(x), Eq.@ and Eq.(|) one obtains 

H(mi, • • • , m-p, 0; x) = ln(x)H(mi, • • • , m^; x) — H(0, mi, • • • , m-p; x) 

— H(mi, 0, m2, • • • , rup-, x) — • • • — H(mi, • • • , r?T,p_i, 0, mp; x) . (21) 

In the case that nip is also zero we can move the last term to the left, divide by two and then 
use again Eq.(|l8|) for all the other terms, thus obtaining an identity which extracts the logarithmic 
singularities due to 2 trailing zeroes. By suitably repeating the procedure as many times as needed, 
we can extract in general all the powers of ln(x) from the generic H-function. A couple of examples, 
if a, 6 are any non-zero indices, are 

H(a,6,0,0;x) = H(0, 0; x)H(a, 6; x) 

- H(0; x) (B{a, 0, b; x) + H(0, a, b; x)^ 

+ H(a, 0, 0, b; x) + H(0, a, 0, b; x) + H(0, 0, a, b; x) , 
H(a,6,0,0,0;x) = H(0, 0, 0; x)H(a, 6; x) 

- H(0, 0; x) ( H(a, 0, b; x) + H(0, a, b; x) ) 
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+ H(0; x) (^H(a, 0, 0, b; x) + H(0, a, 0, b; x) + H(0, 0, a, b; x) 
- (^H(a, 0, 0, 0, b; x) + H(0, a, 0, 0, 6; x) 

+H(0, 0, a, 0, x) + H(0, 0, 0, a, 6; x) ) (22) 



In the same way one can use the product identities, Eq.(|l9|) for extracting the terms singular 
as powers of ln(l — x), or equivalently of i7(l;x) according to Eq.(|l|), around x = 1 from the 
H- functions whose leading (or leftmost) indices are equal to 1. If a = 1 Eq.(p!8|) can be rewritten as 

11(1, mi, • • • , m^; x) = 11(1; x)II(mi, • • • , m^; x) — H(mi, 1, m2 • • • , m^; x) 

— II(mi, m2, 1, • • • m,p] x) — ■ ■ ■ — H(mi, • • • , nip^imp, 1; x) . (23) 

If nil has also the value 1 we can take the second term of the r.h.s. to the left, divide by two and 
obtain an identity to be used when the first 2 indices are both equal to 1 and so on. Let us show 
a couple of examples in the case of two indices a, b not equal to 1: 

H(l,l,a,6;x) = H(l,l;x)H(a,&;x) 

- H(l;x)f H(a,l,6;x) +H(a,6,l;x; 



+ II(a, 1, 1, b; x) + H(a, 1, b, 1; x) + II(o, b, 1, 1; x) 
H(l,l,l,a,6;x) = H(l,l,l;x)H(a,6;x) 

- H(l,l;x)( H(a,l,6;x) +H(o,6,l;x) 



+ H(l; x) (R{a, 1, 1, b; x) + H(a, 1, b, 1; x) + H(a, b, 1, 1; x)^ 
H(a, 1, 1, 1, b; x) + H(a, 1, 1, b, 1; x) 
+H(a, 1, b, 1, 1; x) + H(a, b, 1, 1, 1; x)) ; (24) 



the structure is very much the same as in the equations for extracting the ln(x) singularities related 
to the trailing zeroes. 

It is to be noted that the two procedures - the "extraction" of leading I's and trailing O's - can 
be combined, to give, for instance 

H(l,l,-l,0;x) = ^H(-1;x)H(0;x)h2(1;x) -H(-l,l;x)H(0;x)H(l;x) 

+ H(-l,l,l;x)H(0;x) - ^H(0, -1; x)H(l; x)H(l; x) 

+ H(0, -1, 1; x)H(l; x) - H(0, -1, 1, 1; x) , 

H(l,l,0,0,0;x) = ^H3(0;x)H2(1;x) 

- H(0, 0, 0, 1; x)H(l; x) + H(0, 0, 0, 1, 1; x) 



9 



+ H(0, 0, 1; x)R{0; x)H(l; x) - H(0, 0, 1, 1; x)H(0; x) 
- ^H(0,1;x)H2(0;x)H(1;x) 

+ ^R{0,l,l;x)R\0;x) . 



(25) 



Therefore, one can always express a H- function with leading I's and trailing O's in terms of products 
of powers of H(0; x) and H(l; x), which exhibit the logarithmic singularities in those points, and of 
other "irreducible" H's, i.e. H's whose first index is not 1 and the last index is not and therefore 
is finite at both x = 1 and x = 0. 



We can push further this kind of reduction, by writing all the possible product identities Eq.(|l9| 



and the integration by part identities Eq.(17) and using them for expressing as many as possible H's 
of weight w and "unwanted" indices in terms of products of a "minimal" set of H's of lower weight 
and "accepted" indices. It is to be noted that the number of the H's in the "minimal" set is fixed, 
but their choice is not unique, even if the condition of the extraction of the leading I's and trailing 
O's is imposed. It is easily seen that at weight w the number of relations is nothing but the total 
number of the different products of H's of lower weight and with total weight w. These relations 
are independent when all H-functions of lower weight belong to their respective "minimal sets" . It 
is to be observed, in any case, that the above "reduction" involves only different rearrangements, 
without any modification, of the set of indices which appear in the original H, 

An explicit calculation gives the set sizes of table |. 



Weight 


Full basis 


Irreducible set 


Minimal set 


2 


9 


4 


3 


3 


27 


12 


8 


4 


81 


36 


18 


5 


243 


108 


48 


6 


729 


324 


116 


7 


2187 


972 


312 


8 


6561 


2916 


810 



Table 1: Sizes of the various bases 



The use of the full basis in which each term has only a single H-function gives a unique expression 
in a rather simple way. This is also the preferred representation when higher weights have to be 
built up by successive integration. Expressions can also be given in terms of the irreducible set in 
a relatively easy way. This form is preferred when one has to avoid problems with divergencies. It 
can also be convenient when establishing identities for related arguments. The use of the minimal 
set is particularly convenient for the numerical evaluation of the H-functions, when a large number 
of them has to be evaluated in the same point. It should also be noted that the use of a minimal 
set is relatively easy for the lower weights (at weight 3 it requires only 4 substitutions) while for 
higher weights it will much less straightforward. 



10 



4 Power series expansions 



In general the function }i^{x) does not have a regular Taylor series expansion. This is due to the 
effect that trailing zeroes in the index field may cause powers of ln(x). Hence the proper expansion 
is one in terms of both x and ln(x). Let us first have a look at what happens when there are no 
logarithms. We will use now the other notation for the indices. In that case we have: 



oo j 
X 



■ 1 ^ 

1=1 



H-i(x) = (26) 



i=l 



and assuming that 

^rnix) = E^Snii) (27) 



oo J J 

ax 



1=1 



in which cr = ±1 one can write the relations 

oo j.^ 

i=l * 
oo 

Hl,,7l(x) = V — S'aa,n(«-1) 



i = l 



1=1 i=l 



oo 



lYx^ 



^-l,rh{x) = -V : 5'_^a,n(«-l) 



i=l ' 



= -E^^— ,n^» + E^^n^« (28) 
i=l * i=l 

At this point one could argue what is the better definition of the nested sums. A definition of the 
type 

A Zrf,{i-l) 



n 



-'a,m 

1=1 



E^^^ (29) 



will give only a single term in the expansion and is favored in the mathematical literature, because 
there one is mainly concerned with sums to infinity. For finite values of n however this definition 
has the unelegant aspect that when m has k components that are not zero, the value of Za^mi''^) is 
zero for n < k. We will mostly follow the conventions of ref in which we use the definition: 

SaMn) = E^ (30) 
i=i ^ 



■^Because of the linearity of the problem the presence of more than one term, each with a different Sft would not 
make much of a difference in the following considerations. 
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In this notation one has the property 5^^,(1) = 11^=1 with ai being the sign of rrii. These two 
notations will be referred to as Z-notation and ^-notation respectively. The conversion from one 
notation to the other is not really very complicated if one realizes that Z]}=i = Sj=i Hence 
the 'leading' term has the same index field and the correction terms have fewer indices in which 
some adjacent indices may have been combined. For k nonzero indices there are in total 2'^-! - 1 
correction terms. 

The fact that trailing zeroes in the index field are responsible for powers of ln(x) can be seen 
easily now. Because 

if ^.-In^W = .-.|:(dp_j£>L_ (31, 

we see that once we start with H(Ofc;x), the subsequent integrations due to other indices (the first 
of them not being zero of course, and factors 1/(1 it x) being expanded in x) that come to the left 
of the Ofc will always leave terms with at most k powers of ln(x) and there will be a term with k 
of those powers. Hence the trailing zeroes are responsible for powers of ln(x). Of course the exact 
dependence of ln(x) can be derived much easier by applying Eq.(pl|) repeatedly till all trailing 
zeroes have been removed. This gives an expansion in terms of powers of ln(x) and H-functions 
that are of the type we have just studied and hence can be expanded in x. It is however also 
possible to work one's way through the integrals and the various expansions. This is much more 
work and leads eventually to the same result. Hence we have omitted this derivation. 

If we compare the H-function with the multidimensional polylogarithm in ref we may notice 
that this function can be rewritten into the following expansion: 

oo k j^^j 

= E n-#7 (32) 

with 6o = 1- These functions do not contain powers of ln(6j) and hence they cannot represent all 
H-functions. If we restrict ourselves to H-functions without trailing zeroes one can write the terms 
in the expansion of these H-functions as 

oo k 

E n 7^ (33) 

if we use Z-sums and 



^3 



E n ^7 (34) 



ui>U2>->Uk>l 3=1 ^3 



if we use S'-sums. Hence it is clear that the H-functions without trailing zeroes are special cases of 
the multidimensional polylogarithms with hi = itl/x. For the computation of Feynman diagrams 
we do however need the H-functions with trailing zeroes because of the presence of the logarithms 
(see for instance ref [|lO|] ). 



12 



There is another interesting observation in the expansion. Considering that the expansion of 
an H-function with no traihng zeroes gives terms of the type 



X I 

x=l 

one can introduce another sum by dividing by either 1+x or 1 — x and obtain: 

x=l x=l 

oo 

= {l+x)Y^x\-iyS_^a,A{i) (35) 

x=l 

At times this notation is more convenient. One should however remember that this notation breaks 
down at either 3; = lorat3; = — 1, depending on the particular form used. 

Finally we notice that for x = 1 we have that 

- S.aAoo) (36) 

x=l 

and hence the values of the H-functions in x = 1 are related to the values of the S'-sums in infinity. 
The trailing zeroes do not cause essential problems because when those functions are first written 
in terms of powers of ln(x) these logarithms vanish in x = 1 and we keep only the terms with 
H-functions without trailing zeroes. For the numerical evaluation of these objects one can use the 
algorithms of ref that relate them effectively to combinations of H-functions in x = 1/2 after 
the appropriate conversions. This is particularly interesting for the higher weights because up to 
weight 7, 8 or 9 it is still possible to obtain expressions in terms of a very small number of constants 
(see ref and [^), but beyond these weights this becomes too time consuming while an expansion 
in X = 1/2 is sufficiently fast for nearly all numerical applications, provided that only a limited 
number of them is needed. 



5 The algebra 

The harmonic sums form an algebra |p in which the product of two sums with the same argument 
and having weights wi and W2 respectively can be written as a sum of terms, each with a single 
sum of weight wi +W2- There are two sets of algebraic relations: the relations based on the shuffle 
algebra which hold for all values of the argument, and the relations based on the triangle theorem 
of ref which hold only for values in infinity, provided that not both harmonic sums are divergent. 
For the H-functions we have the general product formula based on Eq.(^). This formula is related 
to the algebra of the harmonic sums, because the harmonic polylogarithms can be expressed in 

^Thus far the only known exact method to do this involves solving simultaneously for all 2 3"'"^ H-functions in 
X = 1. See also a previous footnote. 
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terms of series expansions in which the coefficients are harmonic sums: assume for the moment 
that neither rh nor n have traihng zeroes. In that case we derive: 



in which one of the two powers of 1/(1 — x) has been absorbed in the sum over i. By combining 
the powers of x this formula can be rewritten as 



Note that the inner sum can be done and gives a set of terms that are all single S functions, even 
though the expression may not be very compact. It is called a triangle sum and an algorithm for 



single H- functions. 

For H-functions in x = 1 we have seen that they can be directly expressed in terms of harmonic 
sums in infinity. Therefore the general algebraic rules for those sums that are based on the shuffle 
algebra for harmonic sums can be applied. Hence we see a duality here: the general rules for the 
H-functions correspond to the special triangle rules for the harmonic sums, and the special rules 
for the H-functions in x = 1 correspond to the general shuffle rules for the harmonic sums. 

There is one complicating factor when values in x = 1 are considered. Let us start with assuming 
that the basic divergence H(l; 1) can be used as a symbol. In the case of a 'proper' limit procedure 
such things can be done and after the divergences cancel the finite result should be correct. This is 
called regularization. The general algebraic relations are based on the triangle sums, rather than 
on the shuffle algebra, and the triangle sums are not correct when both objects are divergent. The 
subleading terms will be incorrect. This can be illustrated easily: 




(37) 




(38) 



it is given in one of the appendices of ref |^]. It is also available as a procedure in the language 
of FORM [^. As a result one obtains an expression which can be resummed and gives terms with 




(Hi(x))2 



2Hi,i(x) 





51(00) 




i=l 



25i,i(oo)- 2^2(00)) 




(39) 
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and we see that 



(limHi(x))2 7^ lim(Hi(x))2 . 

X— >1 x^l 



(40) 



The solution to this problem is to be found in S'-space. There it is possible to regularize the infinite 
sums in a consistent way by replacing the sum to infinity by a sum to M with M very large but 
finite, then one can have the divergences cancel and finally take the limit M — > oo. This does not 
correspond to anything one can do in x-space. Because the triangle theorem does not hold for two 
S-sums that are divergent, one cannot apply the regular algebraic relation for H- functions that are 
both divergent in x = 1. Hence the proper algebraic relations at x = 1 have to be derived by means 
of the shuffle algebra which holds for all S'-sums: 



This way is consistent and will allow us to define the Mellin transform properly in one of the next 
sections. It involves the use of values in x = 1. 

Because of the use of different algebraic relations for x 7^ 1 and x = 1, it may happen that 
expressions look rather complicated, but the various algebraic relations between H-functions in 
X = 1 could simplify the expressions considerably. However at the moment there is no known 
systematic method to apply these relations in such a way that one does not have to solve for all 
values in X = 1 first. This way all such objects can be expressed in a minimal independent set of 
objects. Unfortunately there are very many of these objects for a given weig ht w (2 S"'-^) and 
even more relations and hence it is a formidable task to determine all values at x = 1 in terms of a 
minimal set of constants when the weight is large. If the final answer is supposed to be finite one 
can however extract the powers of the basic divergences (they correspond to leading indices that 
are 1) and hence still obtain a finite answer that can be evaluated numerically. The coefficients of 
the divergences can be checked to be zero numerically as well. 

6 Identities between H-functions of related arguments. 

In this section we will look at the identities which can be established for suitable changes of the 
argument. The common feature is that any H-function of weight w and argument x can be expressed 
as an homogeneous expression of the same weight w, involving either H-functions depending on 
a same argument, say t, related to x by the considered change, or constants corresponding to 
H-functions of special constant values of the arguments (typically 1). 

The simplest change of the argument is the change x — > — x. We have seen its effect already in 



Next is the relation between H-functions of x^ and of x. Because 1 + x^ is not a particularly 
interesting object we will have to exclude indices equal to -1 in the H-functions of x^. Restricting 




2S'i,i(oo) - 52(00) 




(41) 



Eq.dD. 
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the indices to only 1 and 0, we can proceed recursively on the weight. For weight 1 we have from 
Eq.i): 

H(0;x2) = 2H(0;x) 

H(l;x2) = H(l;x)-H(-l;x) , (42) 

so that the H's of argument are expressed in terms of H's of argument x, as required. 
For It; > 1, if = 0^, 

H(0^;x2) = 2"'H(0^;x) ; (43) 

otherwise, if = (a, m^_i) for the two cases a = and a = 1 we have, by using the change of 
variable x' = t'^ 



H(0,m^„_i;x^) = / — H(m^_i;x') 

= 2 r^H(m^_i;t'2) (44) 
JO t' 
rx^ fix' 

H(l,m^„_i;a;^) = / -H(m^_i;x') 

Jo I — x' 

The expression of the H(m^_i; t'^) in terms of H's of the same weight and argument t' is supposedly 
known (as we proceed recursively on the weight w); by substituting such expression and then using 
the very definition Eq.(^) all the required j;^ — > x identities are obtained. An example of weight 
w = 2 is 



H(0,l;x2) = 2r^H(l;t'2) 
Jo t 

= 2 r^(H(l;tO-H(-l;t')) 
JO t 



'V 

2H(0, l;x) -2H(0, -l;x) (46) 



and Eq.(15) leads to the well known relation Li2(a;^) = 2Li2(x) + 2Li2(— x). We can observe here 
that a limited set of x'^ ^ x identities could be written only for the Nielsen's polylogarithms 
corresponding to the H„(x) in the notation of Eq.(0), while for the hpl's the set is wider; as an 
example, one can derive for = 3: 



H(l,0,l;x^) = 2(^H(l,0,l;x)-H(-l,0,l;x)-H(l,0,-l;x)+H(-l,0,-l;x)j (47) 

The next transformation of the argument we consider is x — > 1 — x which applies again to a 
smaller set of Nielsen's polylogarithms. Like the previous transformation it is of interest only when 
there are no negative indices (1+x 2 — x is not something we can work with). Proceeding 
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recursively on w, as before, for w = 1 we have 



H(0;l-x) = -H(l;x) 

H(l;l-x) = -H(0;a;). (48) 

The extension to higher weights requires a minimum of care. H(a,m^_i; 1 — x) of weight u; > 1, 
with the first index a equal to or to 1 is the generic function. As discussed in Section 3, if a = 1 
the function can be expressed in terms of a reduced set of functions, where the leading index 1 is 
carried only by H(l; 1 — x), for which Eq.(|4^) holds; therefore, only the case in which the first index 
a is is to be considered. In that case, the change of variable x' = 1 — t' gives 

H(0,m^_i; 1 - x) = / — H(m^_i;x') 

Jo x' 

= / — ;-H(m^_i;x') - / — H(m^_i;x') 

Jo X' Jl-x X' 

= H(0,m^_i;l) - / -H(m^_i;l-t') , (49) 

Jo t — I 



where the constant H(0, rn^_i; 1) is finite (it can be observed here that if the first index is 1 one 
runs into the problem that H(l,mtt,_i; 1) could be divergent). In the general case }i{rhw-i; 1 — t') 
will not be irreducible. We can express it in terms of the H's of an irreducible set of weight w — 1, 
use the supposedly known x = 1 — t identities of weight w — 1 and finally obtain the required weight 
w identity by using the definition Eq.(^). As an example we have at weight 4 

H(0,0,1,1;1 -x) = H(0,0,1,1;1) -H(l;x)H(0, 1,1; 1) +H(1, 1,0,0; x) 
= H(0,0,1,1;1) -H(0,l,l;l)H(l;x) - H(0, 0, 1, 1; x) 

+ ^h2(0; x)h2(1; x) - H(0, 1; x)H(0; x)H(l; x) 

+ H(0,0,l;x)H(l;x) +H(0, l,l;x)H(0;x) . (50) 



A transformation which applies to all the Nielsen's polylogarithms, Eq.(16) is 



x = l/y ; y = l/x ; (51) 

it will be shown that it applies as well to all the H-functions. Before continuing, let us recall that 
the Nielsen's polylogarithms have a (logarithmic) branch point at x = 1, but are otherwise analytic 



for smaller values of x, including all the negative real axis; for studying the transformation Eq.(51) 
it can be therefore convenient to establish the identities for negative values of x, and then continue 
analytically to positive values. The analytic properties of the H-functions are more complicated. 
First of all, if the rightmost index is equal to 0, they have a branch point at x = 0; that is not a 
problem, as we have already seen that we can express any H-function in terms of the functions of 
a reduced set where the trailing index is carried only by powers of H(0; x) = Inx, whose analytic 
properties are well known. If the rightmost index is not and all the indices are in general equal 
to 1 or 0, the H-functions have the same analytic properties as the Nielsen's polylogarithms; but 
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if some of indices are equal to —1, a branch cut at x = —1 appears. Therefore, in the general case 
when indices equal to —1 are also present (and that is the case even of the reduced and minimal 
sets, see Section 3), there is no advantage in considering negative values of x, so that we will start 
from the beginning with an argument equal to x + ie, where x is real and satisfies the constraints 
< X < 1, while e is positive and infinitesimally small; correspondingly, 

y = l/x-ie, (52) 

i.e. the real part of y is also positive, but y > 1, while its infinitesimal imaginary part is negative. 

As in the previous cases, we will proceed by induction on the weight w of the H-functions. At 
w = 1 we have 

H(0;y) = -H(0;x) , 
H(l;y) = H(l;x) +H(0;x) - m , 
H(-l;y) = H(-l;x)-H(0;x) ; (53) 

the constant vr has appeared; it must be given weight 1, so that all the formulas will remain 
homogeneous of the same weight. When continuing the above equations to negative values of x, in 
the interval —1 > x > 0, H{0;x) = ln(x + ie) will develop a positive imaginary part; in particular, 
one has 

H(0;-l)=i7r, (54) 
so that H(l; —1) takes the real value — In 2, as expected. 

For w > 1, rfiw = (a, m^_i), we can proceed by induction, along the following lines 

ry 

B.{m^;y) = / dy' f{a;y')B{mw-i;y') 
Jo 

ry 

dy'f{a;y')B{mw.i;y')+ dy' f{a;y')B{my,^i;y') 
= H(m^;l)+^ ^/ (^a, ^) H (^m^_i; ^) . (55) 

It is to be noted that one can assume that the first index a is different from 1; indeed, as seen in 
Section 3 any H-function of the form H(l,m^_i;?/) can be expressed in terms of a reduced set of 
functions, where the leading index 1 is carried only by powers of H(l; y), whose transformation is 
given by Eq.(^). For a different from 1, H(m^; 1) is a finite constant and the above formulae are 
meaningful. One further finds 

x'2 ( ' x' ) ~^ J x' ' 







-f(-i;-) 

x'2 V ' x' J 



1 + x' 



substituting in the r.h.s. of Eq.(p5[) the identities (of weight w — 1, and therefore known in an 
approach by induction) which express H(m^_i;y' = 1/x') in terms of H(m'^_i;x'), one obtains a 
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combination of terms of the kind 



/ dx' /(a;x')H (x') = H(a, m'^_i; 1) - H(a, m'^„_i; x) , 
and the identities of weight w are estabhshed. As an example, we give the w = 3 identity 
b(o,-1,1;- -ie] = -H(0, -1, 1; x) + 2H(0, -1, 1; 1) 



+ 2H(0, 0, -1; x) - 2H(0, 0, -1; 1) + H(0, 0, 1; x) - H(0, 0, 1; 1) 
H(0, -1; x) + H(0, -1; 1) + H(0, 1; 1)^ H(0; x) + ^H='(0; x) 

-ivr (^H(-l;l)H(0;x) + ^H2(0;x) -H(0,-l;x) +ii'(0,-l;l) ) . (56) 

Another important set of identities, which is however vahd for any set of indices and has 
no counterpart within the Nielsen's polylogarithms, applies to arguments x and t related by the 
transformation 

x=i^. (57) 

whose inverse is again 

t = ^ . (58) 

1 + x ^ ' 

Even in that case, it turns out that any H-function of weight w and argument x can be expressed 

as a homogeneous expression of weight w, involving H- functions of argument t, related to x by 

Eg. (|57|) , as well as constants corresponding to H-functions of argument 1. The proof is, again, by 

induction on the weight. If w = 1, from the very definition Eq.(||) one immediately finds 

H(0;x) = -H(l;t)-H(-l;t) , 
H(l;x) = -H(0;i) -H(-l;l) +H(-l;t) , 
H(-l;x) = H(-l;t) -H(-l;l) . (59) 

For w > 1 and = 0^ the result is trivially true, as can be verified by inspection; the same is 
true also for = Iw and rhw = —Iw In the more general case, write rhw = (a,m^_i); where the 
index a takes the values 0, 1, —1. As discussed in Section 3, and already recalled for the x — > 1 — x 
identities, if a = 1 the function can be expressed in terms of a reduced set of functions, where the 
leading index 1 is carried only by H(l; x), for which Eq.(|59|) holds. In the other two cases a = 0, — 1 
the change of variable 

1-t' 

X — 



l + t' 



gives 



dx 

H(0,mu,_i;x) = / — H(m^_i;x' 
Jo x' 
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/■I dx' 

H(0,m^„_i; 1) - / — H(m^_i;x') 

J X X 

H(0, 1) - dt' + ^) H 



H(-l,m^_i;x) = / — — -H(m^_i;x') 
JO l + x' 

= R{-l,m^_,;l) - dt'j^u(m^^v,]-^^ • (60) 



At this point, one can substitute the relations already found to be valid at weight w — 1, for 
expressing the functions H (m^_i; (1 — t')/{l + t')) in terms of H's of weight w — 1 and argument 
i', and then perform the last integration in t' according to the definition Eq.(|6|). 

As an example, we give the following w = 3 identity 

H(-l,-l,l;a;) = -H(0, -1, -1; i) + H(-l, -1, 1; 1) 

+ H(0, -1; t)H(-l; t) + Jh3(-1; t) - ^H2(-1; t)H(0; t) 

-iH(-l;l)H2(-l;f)-H(-l,l;l)H(-l;t) . (61) 

7 Identities between H's and related functions. 

Let us introduce a related set of functions G{fhw;x), where has almost the same meaning as 
for the H's, but the first index is always equal to 1, i.e. m = (1, m^_i), through the definitions 

/■I df 

G(l;x) = -/ --- (62) 
Jo t - 1/x 

for w = 1 and 

/•I dt 

G(l,m^_i;x) = - / - — --—R{myj-i;t) (63) 
Jo t - 1/x 

for w > 1. 

The G(m^;x) are nothing but homogeneous combination of H- functions of weight w. As by 
now usual, we will show it proceeding by induction on w. For w = 1, hy performing explicitly the 



elementary integration we obtain from Eq.(62) 



G(l;x) =H(l;a;) . (64) 

Next, assume that the identities are established for w; put fh = (a, r7i^_i), and consider the 
functions of weight w + 1 given by 



ri cLt 

G(l,a,m^_i;a;) = - / - — -— H(a, m^_i; t) 
Jo t-l/x 



(65) 
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One can differentiate with respect to x, then integrate by parts in t, using of course Eq.(|63|) when 
relevant; considering for instance the case a = —1 one obtains 

d 

— G(l,-l,m^_i;x) = [f{-l,x)-f(0,x)]G{l,mw-i;x) 

+ [/(-l,x) + /(l,x)]H(-l,m^_i;l) . (66) 

Similarly, one has 

d 

— G(l,0,m^_i;x) = -/(O, x)G(l, m^_i; x) 

+ /(-l,x)H(0,m^_i;l) 

^G(l,l,m^_i;x) = [/(0,x) + /(l,x)]G(l,m^_i;x) . (67) 

One can substitute the already obtained identities expressing G(l,m^_i;x) in terms of H's of 
weight w and then integrate in x between and x by using the very definition Eq.(^) (according 
to Eq.(|6^) the G-functions vanish at x = 0). The required identities of weight w + 1 are then 
established. As an example, we give one of the identities of weight w = 4: 

G(l,0,-l,l;x) = -H(0,-l,0,l;x) -H(0,-l,l,l;x) 

+ H(0,0,0,l;x) +H(0, 0,1,1; x) 

- H(-l, 1; 1)H(0, -1; x) - H(-l, 1; 1)H(0, 1; x) 

+ H(0,-l,l;l)H(l;x) . (68) 

In the same way one can work out the similar identities existing for several related classes of 
functions such as, for instance, 

-R{a;t)R{b;xt) 



or 



t - 1/x 

[\tf{a,t)i{{d;t)i{{b-xt) . (69) 
Jo 

8 Special values of the H's and their numerical evaluation. 

It is known that the Nielsen's polylogarithms for the special values of the arguments equal to 
+1,-1 and 1/2 can be expressed in terms of a few mathematical constants, typically Riemann 
^-functions of integer arguments; the representations which they provide for those constants as 
definite integrals can be manipulated by means of integration by parts, changes of variables and 
the like providing the analytic values of a number of definite integrals of special interest. The same 
applies, and in much more systematic way, to the H- functions, thanks to the greater and wider sets 
of identities which they satisfy. 
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In the case of the H's, it is not necessary to consider as independent the values corresponding to 
the argument equal to —1; indeed, one can always express any H- function in terms of the reduced 
set of functions in which trailing indices equal to are missing, so that by using Eq.(|l5|) one can 
replace a value at x = — 1 with the value at x = 1 of a related function. In analogy with the Nielsen's 
polylogarithms case, it is convenient to consider also the values at x = 1/2 of the functions whose 
indices are equal to or 1 (i.e. when the index —1 is missing). 

More specifically, one can consider: 

• the x^ — > X identities, Eq.s(p2[-^), for x = 1 ; 



the 1 — X ^ X identities, Eq.s([lS|-|50[). They can be used at x = 1/2, providing a first set 
of identities for the values at x = 1/2, but also at x = —1; in the second case, one gets 
values at X = 2, which are converted into values at x = 1/2 by using the x — > 1/x identities, 



Eq.s(51-56), as well as values at x = —1, which are converted into values at x = 1 by Eq.(15); 



the just recalled x — > 1/x identities, Eq.s(51-5(:), at x = 1 and x = —1, followed by the usual 
conversion to x = 1 through Eq.(^); 

the X — > (1 — t)/(l + t) identities, Eq.s(^ |6l| ), at x = corresponding to t = 1 (they are 
automatically satisfied, by construction, at x = l,t = 0); 

one more set of identities is obtained by writing the identities between G-functions and H- 
functions, discussed in Section 8, at the special value x = —1, by using the relation, which 
follows from the definition Eq.(|63[) 

G(l,m;-1) = H(-l,m;l) 

and converting once more the values at x = — 1 of the H's into values at x = 1 by means of 
Eq-dll). 



The set of relations obtained in that way is highly redundant; it has been checked explicitly that 
they generate the table of the u) = 4 definite integrals given in Appendix B of the second reference 
of [^. It has not yet been investigated whether they are sufficient, by themselves, to generate also 
the tables of higher weights obtained in [^] . 

Another powerful method to obtain the values at x = 1/2 when there are no negative indices 
is by considering the transformation x + z), which corresponds to a suitable combination 

of the transformations x ^ 1/x and x (1 — x). Using the same techniques as in the section 
on related arguments, all these objects are directly expressed in terms of H- functions in x = 1. 
Such expressions can then be used in reverse to obtain the numerical values of the 'independent 
constants' that occur in the expressions for the H- functions at x = 1. As an example we have 

^^^G) = ^C3-^C2ln(2) + lln3(2) (70) 
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which is of course weh known. We have also 

Both relations provide a power series for the evaluation of Cs- The method gives also an expression 
of sq = S-5-i{oo) in terms of H5^i(l/2), 'H.q{1/2) and combinations of constants of a lower weight. 
Similar dependencies can be derived for the higher weight constants. 

Let us finish with a few remarks on the numerical evaluation of the H's for arbitrary values of 
X. According to the discussion of Section 3, it is sufficient to restrict ourselves to the H's either 
of the reduced set or of a minimal set, as all the others can be obtained from them as suitable 
combinations. The H's of such a set have no trailing indices equal to 0, so that they can be expanded 
in series of x around x = 0. For small values of x the series will be rapidly convergent, but the 
convergence will slow down approaching the cuts at x = ±1. But for x approaching 1 we can use 
the transformation Eg. (|57|) , so that the corresponding t = {1 — x)/{l + x) will fall in the region 
near and the expansion in t will be rapidly converging. 

More exactly, the equation 

1-r 

r = (72) 

has the two solutions r = —1 — and r = — 1 + ^/2. Therefore, we can use the expansion 
around a; = in the interval — (\/2 — 1) < x < \/2 — 1, where |x| < \/2 — 1 < 1/2, switching for 
V2-1 < X < -v/2 + 1 to t = (1 — x)/(l + x), which corresponds to \t\ < — 1 < 1/2. For greater 
values of x, one can use the x — > 1/x identities. For large negative values of x, i.e. x < 1 — -v/2, 
one can flip the sign of x with Eq.(|l5|) and then proceed as above. 

In practice the transformation of Eq.(57) can lead to a large number of functions to be evaluated 



and hence it may be more profitable to apply this transformation only for values of x that are much 
closer to one. If, on the other hand, nearly all H-functions of a given weight have to be evaluated 
for some value of x one can use the turnover value of \/2 — 1 in a rather profitable way. 

The values in x = 1 require some extra attention. These are actually needed rather frequently 
and hence there exists some literature on them. From Eq.(p8|) it should be clear that an H- function 
in X = 1 can be expressed in terms of either S-sums or Z-sums in infinity. Hence much information 
can be found in |||] and the papers they refer to. Ref. gives a different method to evaluate 
these sums. Recently this method has been used by one of us (J.V.) to obtain all such sums up 
to weight 9 (see also footnote 2). For only nonnegative indices results have been obtained up to 
weight 11 |]lT|. When the first index of the H-function (or the 5-sum) is one, the value in x = 1 (or 
the sum in infinity) will be divergent. Yet we have to consider these objects. As mentioned in the 
section on the algebra this can be done consistently only in terms of the sums. Hence the safest 
method is to rewrite the H-functions in x = 1 immediately in terms of either 5-sums or Z-sums. 
In the case that the weights are low enough, these can then be rewritten in terms of a limited set 
of 'fundamental constants'. 
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9 Mellin transforms 



At times one may need the Mellin transform of the Harmonic polylogarithms. In ref a method 
is given to evaluate such transforms for a class of functions which is more or less the class of H- 
functions. There is however one complication with Mellin transforms. Divergencies at rr = 1 must 
be extracted. This is because the Mellin transform is defined by 

M{f{x),N) = f\xx^f{x) 



(1— JO 1 — x 

^ ^ -/(!)) in^d-.) 

in which the function / is supposed to be finite for x = 1 when the factor 1/(1 — is present. 
Hence we have to pay attention to the powers of ln(l— x). They can be isolated with Eq.(|2^). After 
this extraction the remaining H- functions are finite in x = 1. 

At this point we can attack the Mellin transforms. It is easy to obtain the lowest weight results: 

1 1 

dxx''R{0;x) = 



n + 1)2 

\xx-Bil;x) = ^li^ 
n + 1 

/'d.."H(-l;.) = Hr^.l(^ + M?l,i + H,..) p4) 

JO n + 1 n + 1 

in which we have used that H(— 1;1) = — 5_i(oo) = ln(2). The higher weight results can be 
obtained by recursion. Like in ref this is done by partial integration. We also exchange the sums 
immediately after each step so that we may do one of them immediately. The result is: 

/ dx (7'x'Ho,^(x) ^^^^~^^k = O-Ho^m(l) (^.^(fc+l),?!^^) - 
^ i=n ^ ' 

^xf:aVH^(.)M±l) (75) 



S^{i+l) 



1 oo 

/ cix^x*H^(x)(^o-S'^(fc+i)_p'(i+l) 

^ i=n 



-cr ..^ - crS^(^k+i),p{n) j (76) 



(i + 1)^+1 



dx J2 cr'x'B_i^rnix) ^(f^~^^k = o-H_i,^(l) (5'^(fe+i),p'(oo) - 5'^(fc+i),p-(n; 
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„1 oo 
■ „ 



The variable a is either 1 or —1. This leaves only the evaluation of the H- functions in x = 1. These 
values do not have to be finite. Only the H- functions that are used in the subtraction in Eq.(|73[) 
are finite. This causes no problems provided the divergencies are regularized in the representation 
in terms of S'-sums as explained before. 

As an example we show here a nontrivial Mellin transform: 

^(%^r^'^) = 5l,-2,-l,2(iV)-25l,_2,-3(A^) + 2Si,5(iV) 

-i5i,_2(iV) (^(2 ln(2) + C3^ + Si4N) (^^C2 ln(2) - Cs) 
+5i,i(iV) (aLu ln^(2) - C2 ln2(2) - 




.27 6 

+Si{N) (^IC2C3 - f C5) - ^C2 ln^2) + ^C2C3 ln(2) 

-C2Li4(^) + ^C|ln^(2) + ^Cl - ^Cl + ^5-5,-1(00) (78) 

The sum in the last term is irreducible. 

In the case that the weight of the terms too large (currently larger than 9) it becomes rather 
hard to obtain the values for the H-functions in x = 1 or alternatively for the S-sums in infinity. 
Because the algebra for the H-functions in x = 1 is different from the algebra for the H-functions 
for general values of x there may be large numbers of H-functions left that each are divergent at 
X = 1. The reason is that some algebraic work is done first with the general algebraic rules and 
has to be 'undone' with the rules for x = 1. The relations that make the divergences cancel may 
not be easy to find. One can still obtain numerical results however. 

If one is faced with higher weights one may proceed as follows. The H-functions in x = 1 are 
first expressed in terms of S-sums in infinity. Then the shuffle algebra for the S'-sums is used to 
extract the divergencies in a way that is similar to how this is done for the powers of ln(l — x) 
for the H-functions. Because the divergences have to cancel each other, all divergent terms should 
disappear, even though we may not have the algebraic methods to prove this for the case at hand. 
The remaining finite expression can in principle be evaluated numerically. 

Inverse Mellin transforms are now relatively easy. As pointed out in ref each S-sum has 
a single most complicated original function in terms of H-functions in which we can define 'most 
complicated' by function with largest weight or in the case of identical weights the largest number 
of nonzero indices. And actually one can obtain the relation between the S-sum of which one needs 
the inverse Mellin transform and this most complicated H-function from the recursion relations in 
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Eq.([75|). Hence the algorithm is clear: 



• Locate the most complicated S'-sum(s). 

• Construct the corresponding H-function(s) in x-space. 

• Add it and subtract it. 

• Make the Mellin transform of the subtracted version. This will cancel the original S'-sum. 

• Repeat the above steps until there are no more S'-sums remaining. 

• Multiply the remaining constant terms by 5{l — x). 

This algorithm will properly terminate. It has only one problem: Some Mellin transforms have a 
factor (—1)^ and some don't. What if we take an S-sum which should have a factor (—1)^ but we 
omit it? Here we have to realize that the inverse Mellin transform is to be constructed from either 
all even or from all odd moments only. Hence we have to specify whether is even or odd. This 
will give a value to (—1)^- Hence the only thing that remains is to give the relation between an 
S'-sum and the most complicated H- function that contributes to it. 

• If the number of negative indices is odd, there will be a factor l/(l+x), otherwise there will 
be a factor 1/(1 — x)+. 

• Next copy the index field to the H-function. 

• Working from the rightmost index to the left, each index will get a sign that is the combination 
of its old sign and the signs of all indices to the left of it. 

• There will be an additional overall sign on the term which is the sign of the last index. 

• There will be an additional overall sign on the term which is (— l)"'~i~'^ in which w is the 
weight of the S'-sum and d its depth (which is the number of nonzero indices) . 

• Each negative index in the current configuration will give a minus sign to the term. 

We will give two examples of weight 7 functions. First and example that involves subtractions with 
ln2(l — x) in the Mellin transform: 




(79) 
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In this case there is no difference for even values of N and for odd values of A'^. However the next 
example is different. For even values of A'^ we have 

-S-1,1 -2,1,2 (-/V) ^ ^-H_i _i,2,i,o(a;) + 2H_i _i,o,o,o,o(a:;) - 2Ho,o,o,o,o,o(a;) 

+^H_i,_i(x)C|-Ho,o(x)C| 

1 67 
+H_i(x)(-C2C3 + C|ln(2) + -C5) 

1 67 
+H„(x){--C2C3 - 2Cl ln(2) - -Cs) 

-l<.(sH2) - |C| ln^(2) - |IC| + id + |C. ln(2) - 5..) 

+^((H,(x) + 2H„W)(l6C3-|« 

61 n 35 ^2 93^ , 3 \ 

-560^^ + T28^3 + 64^^^^(')-r^j 

+6{l-x) (-^C2C3 ln^(2) - UC2C5 + Y^C2 ln^(2) - CsLisC^) 

__93_ 2a —C'^\n^(2) - -^^^111(2) - ^^^C'^ln(2) 

140^^^^ 12^^ ^ ^ 280^^ ^ ^ 896 ^ ^ 

197, , o,^, 37215, 19, 10 29 \ 

Csln (2)H CtH ln(2)cr6 cryaH a-j b (80) 

S5 V ; -r 35g4 "-^ ^ 28 ^ 7 ^ 14 '^''^ y ^ ^ 

in which 

(76 = 'S'_5_i(oo) 
0"7,o = 'S'_5,i,i(oo) 

(77,6 = ^5,-1,-1(00) (81) 

In the case of odd N the terms with 1/(1 + x) change sign. As one can see these formulae can 
become rather involved, even though the number of terms is rather small compared to the number 
of functions that exist in x-space for this weight. 

In the case that sums of a higher weight are considered one may not be able to substitute the 
values of the H- functions at x = 1. The same considerations as for the Mellin transforms can be 
used to obtain an answer that can at least be evaluated numerically. In general the formulae will 
of course be much lengthier. 
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